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Abstract 

We derive the (matrix-valued) Feynman rules of mass perturbation theory of the 
massive Schwinger model for non-zero vacuum angle 9. Further, we discuss the properties 
of the three-boson bound state and compute - by a partial resummation of the mass 
perturbation series - its mass and its partial decay widths. 
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1 Introduction 



The massive Schwinger model is two-dimensional QED with one massive fermion. In 
this model some nontrivial field-theoretic features occur. E.g. instanton-like gauge field 
configurations are present in the model and, therefore, a ^-vacuum has to be introduced as 
a new, physical vacuum (0, §]). Further, confinement is realized in this model and there 
are no fermions in the physical spectrum (0, m^)- The fundamental physical particle of 
the theory is a massive, scalar boson (the Schwinger boson with mass fi). This boson is 
interacting and may form bound states (0, |^, |^) and undergo scattering (0). A proper 
way to discuss the model is mass perturbation theory (i.e. perturbation in the fermion 
mass term), as it preserves the nontrivial structure of the model (like e.g. 6'- vacuum) 
and is performed in terms of physical (bosonic) variables only ([|], |], |10], [ll|]). This mass 
perturbation theory uses the known exact solution of the massless Schwinger model as a 
starting point (in the massless model the Schwinger boson is a free scalar particle with 
mass /Uo = 
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In this paper we first derive the Feynman rules of the mass perturbation theory 
for general vacuum angle 6 (which slightly complicates them as they acquire a matrix 
structure). We then discuss the three-boson bound state and show how, by a resummation 
of the perturbation series, we may compute the mass M3 of the three-boson bound state 
and its partial decay widths for the decays into two Schwinger bosons and into one 
Schwinger boson and one two-boson bound state. 

The mass perturbation theory that is used throughout this paper is based on the 
Euclidean path integral formalism and, therefore, we have to take into account our specific 
Euchdean conventions (see e.g. |T7||). 



2 Mass perturbation theory 

First let us shortly review the mass perturbation theory. By simply expanding the mass 
term, the vacuum functional and VEVs of the massive model may be traced back to 
space-time integrations of VEVs of the massless model. E.g. the vacuum functional is 



Z(m, 



00 ^„ n 



n=0 i=l 



E n / dx,^!{x,)^>{xi) exp / dx ^(z^ - e4)^ - \f^,F^^ 



{k ...instanton number). Therefore, one needs scalar VEVs {S{xi) . . . S^Xn))^ of the 
massless model, where S = S± = |^(l±75)^f. Chiral VEVs (5^, . . . ShJo, Hi = ±, 
are especially easily computed, as only a definite instanton sector contributes (see e.g. 
0,0,0), 



^2 



Xj) 



(2) 
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k = (Tj = — ri- 

i=l 

where cxj = ±1 for Hi = ±, D^^ is the massive scalar propagator of the Schwinger boson 
= — ) and S is the fermion condensate of the massless model. 

The Schwinger boson $ is related to the vector current, = ■^e^^d''^, and, there- 
fore, the S and $ VEVs, which we need for the perturbative calculation of massive VEVs, 
are related to the vector and scalar current VEVs of the massless model. Explicitly the 
S and $ VEVs may be computed from the generating functional (which is at the same 
time a VEV for n chiral currents) 



■ exp 



see 



/ o!|/i%2A(2/i)D^„(yi - y2)A(2/2) + li^iZ^-T f dyX{y)D,^{y - xi)] (3) 

1=1 

20| for an explicit computation), where A is the external source for the Schwinger 
boson Observe the {—Y' in the last term of the exponent. As a consequence, whenever 
an odd number of external $ lines meets at a point Xi, the corresponding S-{xi) acquires 
a — , i.e. instead of a 5" = 5"+ + S*- vertex there is a P = 5"+ — S"- vertex. 

From equs. (2), (3) one finds that exponentials exp ±47rD^(,(xj — Xj) are running from 
any vertex Xi to any other vertex xj; however, in order to get an IR finite perturbation 
theory, one has to expand these exponentials into the functions 

E±(x) = e±^"^''o(-) (4) 

(or their Fourier transforms E±{p) for momentum space Feynman rules). This expansion 
is analogous to the cluster expansion of statistical physics. 

The mass perturbation theory may be summarized by the following Feynman rules 



Dno(Xi-Xj) ••• 

E+(.) (Xj -Xj) • • • vAAAAA/ 

mE 



go 



X 



Fig. 1 
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where mS is the bare couphng constant of mass perturbation theory and g the renor- 
mahzed one (including all tadpole-like corrections, see e.g. 0). Strictly speaking, we 
added the 5"+ and contributions of each vertex, therefore the rules, as they stand, are 
true only for the special case 9 = 0. The general 9 case we will discuss in a subsequent 
section, where we will find that the graphical notation of Fig. 1 may be used, but the 
meaning of the individual graphs is slightly changed. 

Using this mass perturbation theory the following features could be shown to hold in 

1. The mass corrections for the Schwinger boson (computable from ($(xi)$(x2))m; 
[pip also occur as corrections to the internal boson lines of the exponentials E±{x). 
Therefore, the bare Schwinger mass /iq may be replaced by the renormalized one, 
fi, in all internal propagators. 

2. In the perturbation series of Fig. 1 there occur, among other terms, strings of E±, 
E±(xi — X2) ■ E±{x2 — X3) . . . When the two-boson part ^^-D'j^{x) is separated in 
the expanded exponential E±{x), (4), these two-boson blobs may be resummed in 
momentum space. 




+ 




+ 




+ ... 
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via a geometric series formula and give the leading contribution to the two-boson 
bound state mass pole (with mass M2). Higher n-boson blobs give only higher order 
corrections to the mass, because the two-boson blob is near its threshold singularity 
at M2 (see §). 

Explicitly one finds the mass pole {s = —p^ > 0, remember our Euclidean conven- 
tions) 

1 



1-i^mScos^-^ ^ 



: arctan 



(5) 



--1 



-1 



with the pole mass 



m2 = V-A2 , A2 
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and residue at the pole 



2A, 



Svr^mE cos 9 



mT, cos 9 jj? 
Analogously, we want to show in this paper that 



(6) 



(7) 
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1. there is a resummation that leads to a //-M2-blob 




Fig. 3 

(this is quite obvious from the Feynman rules, Fig. 1 and from the fact that the 
propagators E± may connect all vertices) ; 

2. two- and three- boson blobs and the above-mentioned /x-M2-blob (Fig. 3) occur 
in such combinations that they may again be resummed via the geometric series 
formula to result in a propagator 



1 — mTi cos 9 f{p) 
and f{p) is the sum of the three above-mentioned blobs, 




Fig. 4 

where the three-boson blob will give rise to the three-boson bound-state mass pole 
M3, whereas the 2fi- and /x-M2-blobs have imaginary parts that will give the partial 
decay widths for the decays of the three-boson bound state M3 into 2/^ and 11 + M2, 
respectively. All other types of blobs that we ignored are unimportant for these 
features. 

Actually, things are a little bit more complicated. As mentioned above, up to now 
the Feynman rules of Fig. 1 are true only for the special case ^ = 0. But this restriction 
is too strong. E.g. the decay M3 — > 2// is not possible in that case because of parity 
conservation. So let us discuss the general 9 Feynman rules as a next step. 



3 Feynman rules for general 6 

In the mass perturbation theory each vertex is mS = m(S'+ -|- S-). For general 9 
the two chiral components S+ and S- couple differently, with (renormalized) coupling 
ge — m^e*^ -I- o(m^) for S+ and — m^e~'^ -|- o(m^) for S^. As a, consequence, the 
Feynman rules of Fig. 1 acquire a matrix structure. More precisely, the wavy line (the 
E± propagator) turns into a matrix. 

Sip) = ( ^Ap) I-(^) ^ (9) 
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whereas the vertex g turns into an n-th rank tensor Q when n wavy hnes £ meet at 
this vertex. Only two components of this tensor are nonzero, namely ^4.+...+ = ge and 
Q = gl- E.g. when two lines meet at one vertex, then Q is the matrix 

External bosons $ that meet at a vertex are represented by P (S), when an odd (even) 
number of bosons meets at the vertex, where 

Among other Feynman graphs, there exist the following ones (that may be resummed), 
where we amputate external bosons. 



Fig. 5 

or 



with the solution 



gn{p) := g + g£{p)g + g£{p)g£{p)g + ... 

= g{i + £ip)gu{p)) (12) 



''^''>-N(p)[ goE4p) l-geE^{p)) ^^^> 



where 

N{p) = det(l - g£{p)) = l-{ge + 9*e)E+{p) + ge9*e{El{p) - El{p)). (14) 

So we succeeded in inverting the propagator n(p), analogously to our above discussion. 

For the determination of the mass pole and decay widths only the denominator N{p) 
of (13) will be important. In the denominator N{p) we find all the n-boson propagators 

(47^)" — 

dn{p) := ^-^D-{p). (15) 

But we want to include the M2 two-boson loop 

, , r dPq Svr^mScos^^ An 

'''^^^ J (27r)V'(9' + M|) (p-g)2 + ^2' y ^) 

too. In principle, its existence is obvious from the Feynman rules of Fig. 1. But how does 
it fit into our matrix scheme? The answer is that it is parity odd and may be treated 
just like any odd n-boson propagator (e.g. like ^3), as may be checked by expanding 
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the perturbation series into the individual n-boson propagators and including the /X-M2 
two-boson propagator term by term. So we just include it into E± with a ± sign as a 
higher order contribution that wc need. 

Inserting the n-boson propagators <!„, we may rewrite the denominator N{p) like 

N{p) c^l-{g0 + g;){di + d2 + + d3 + ...)+ 

^9eg*e{di{d2 + 4 + • • •) + «^i,i(c?2 + c?4 + • • •) + c?3(c?2 + (i4 + •••) + •• •) (17) 

where we included di i, as stated above. 

Now suppose we are at the three-boson bound-state mass M3. Then the real part 
of (17) vanishes by definition, and the leading order contribution (of order ^^^^) stems 

from the d^{s — —p^) function that is near its threshold singularity, 

{ge + 9*e)d3{s = M|) = 1 + o{m^). (18) 

In addition, all functions di that are above their respective mass thresholds will have 
imaginary parts at Af|, and result thereby in partial decay widths. This is obviously so 
for d2, and it is also true for di^i provided M3 > /i + M2 (actually this is the case as we 
will see in the computations). Therefore, at M3 (17) may be rewritten like 

-i{9e + 9*e) (im rf2(M|) + Im di,i(M|)) + Aige9*gds{Ml)lm d2{M^) + o(m') 

= -imEcose(lmd2(M|) +Imdi,i(M|)) + i^^Imd2(M|) + o(m^) (19) 

COS u 

and we find that parity forbidden decay widths acquire a factor (cos^ — ^^), whereas 
parity allowed decays acquire a cos 9 factor. 

So we have, indeed, succeeded in rewriting the mass pole equation, as indicated in 
Fig. 4, and may start the actual computations. 



4 Bound state mass and decay width computations 

For the computation of the three-boson bound-state mass M3 we need the three-boson 
propagator ^3 and find, in lowest order (see (18)) 

1 = ^mE cos ■ QAir^Dlip) (20) 
or, after a rescaling p — > ^ to dimensionless momenta 

cosODlip). (21) 



6 //2 



D^{p) is given by the graph (where we introduce positive squared momentum s = —p^ > 
0) 
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p\ rx r 

-2 dx dy 
Jo Jo J 



Fig. 6 

d^qid^q2 1 



(27r)4 (p + g, + g2)2 + lgf + + l 
d^qid^q2 1 



(27r)^ 



r dx 



?i + 1 + (?2 - Qi)x + {{p + qi + q2Y - ql)y 

dy ^ 

s{xy — x'^y — + xy"^) — x + x"^ — xy + y'^ 



r dx 
J 87r2(l - s(l 

1 (ia; 



dz 



87r2(l-s(l-x))T(s,x) 



x)) Jo + T'^{s, x) 

1 x 
arctan 



2T{s,xy 



T'{s,x)^ 



x^ — sx^il — x) + 4x(l — x) 



(22) 



(23) 



4(s(l-x)-l) 

where, as usual, wc introduced Feynman parameter integrals and performed the momen- 
tum integrations. Further, the first Feynman parameter integral could be done analyti 
cally. The numerator of has a double zero at s = 9: 



9x(x--)2. 



(24) 



This double zero is in the integration range of x and is precisely the threshold singularity. 
Setting 

s = 9 - A3 (25) 

in the numerator of in the factor and s = 9 everywhere else, where it is safe, one 
arrives at: 



27r2 L 



1 d.i 



arctaii 



Un'^Jo ^|9x - 8| ^(x - |)2x + ^x2(l - x) 
The mass-pole equation reads 

647r^ 

1 = ——mJ: cos 91 {A3) 



-■■ AA3). 



(26) 



(27) 



and must be evaluated numerically. It gives rise to an extremely tiny mass correction A3. 
For sufficiently small m it is very well saturated by 



A3(mEcos^) ~ 6.993 exp(- 



0.263 



mE cos 9 ' 



(28) 
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and is therefore smaller than polynomial in m. (I checked the numerical formula (28) for 
30 < n^scose < 1000' corresponding to lO^^ < A3 < 10-^°°, but I am convinced that 
it remains true for even larger however, there the numerical integration is quite 

difficult because of the pole in (26).) 

We conclude that the three-boson bound state mass is nearly entirely given by three 
times the Schwinger boson mass (we change back to dimensionfuU quantities now), 



M| = 9u2-A3 , A3 = 6.993u2exp(-0.263— f -) (29) 

mh cos 

or, differently stated, that the binding of three bosons is extremely weak. 

Therefore it holds that M3 > /i + M2, as indicated above, and, consequently, a decay 
of M3 into fi + M2 is possible. This has the consequence that the three-boson bound state 
is unstable even for 6' = 0, contrary to some earlier conjecture (0). 

So let us turn to the decay width computation. For this purpose we need the first 
Taylor coefficient of the denominator 1 — mTj cos 9 d^{s) (see (14)) around the mass pole. 
Observe that because of formulae (25), (28) mS cos9d^{s) may be written, in the vicinity 
of s = M|, hke 

V I \ ^^cosO 6.993/^2 

mL cos 6*03(5) ~ m — . (30) 

0.263 9/i^ — s 

Therefore, we find the Taylor coefficient 

mS cos 6 

C3 = 31 

0.263A3 ^ ^ 

Generally, a decay width may be inferred from the imaginary part of a propagator, 

^/ N const. , , 

G{p) — — -, (32) 



M2 - iMT ' 



where F is the decay width. 

In our case we have for equ. (19), 



Nip) cs{s - Mi) - imS(cos6' - ^)Imd2(M|) - imScos6'Imc/i,i(M|; 

const. 



[cose- -^)lmd2{Ml) + cos^Im(ii,i(M| 



(33) 



Next we need the imaginary parts lmd2, Imiii^i. Both of them stem from a two-boson 
blob, so let us write down the general result (which is standard) (s = —p"^) 



lmiD~DM.)is) = I^/(^ ^.;Mf(p-g)^VM| 

^ (34) 



2w{s,MlMi 
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w{x,y, z) = {x + y +z — 2xy — 2xz — 2yz) 'K (35) 
Therefore we can write for (34) (see (7) and (16) for the normahzation factors of d2 and 

const. 

Q _ A/f2 - Vin^rprx^ 1 'I 47r2 - (mS cos 0)2 l&y^ ^ > 

^ ^^^3 ''c^y'-^^^ cose )w{Ml,fi'2,ii2) ^ C3 fi'2w{Mlfi'2,M:j) 

There seems to be something wrong with the sign of the parity forbidden partial decay 
width TM:i^2n (the d2 term). Actually the sign is o.k. and the problem is a remnant of 
the Euclidean conventions that are implicit in the whole computation (see e.g. |T^, pi]|). 



In these conventions 6 is imaginary and therefore (cos^ — ■^^) > 0. Of course, this is not 
a reasonable convention for a final result. When performing the whole computation in 
Minkowski space and for real 6, roughly speaking, the roles of £'+ and are exchanged 
in (14). This gives an additional relative sign between parity even and odd n-boson 
propagators and, therefore, changes the factor of d2 to {-^^ — cos 9), which is > for real 

e. 

With this remark in mind, and expressing the final results for real 6, we find, by using 
the approximations 

w{Ml, Ml, I?) ~ ^(V, /i^) = 2^3/1 ^/a^ + o(m2), (38) 
the following results: 

.47r2A3, 1 



0.263 



9v^/i cos2^ 



3.608/i(^— - 1) exp(-0.929 ^ ^ (39) 
cos^ o m cos tl 



and 



^M3^fi+M2 — 0.263 

3v3/i 

~ 43.9;uexp(-0.929 — ^) (40) 

mcosp 

where we inserted the numerical value S = ^ = 0.283/i. 

The ratio of the two partial decay widths does not depend on the approximations 
that were used for the M3 computation. 



(41) 



Contrary to a naive expectation it holds that VM■i-*^J.+M2 > ^Ma~^2^i, although /x + M2 ~ 
M3. This is not surprizing in two dimensions, because there the phase space "volume'' 
does not grow with increasing momentum. 
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5 Summary 



So we have succeded in computing the three-boson bound state mass and partial decay 
widths via (resummed) mass perturbation theory. The computations may be generahzed 
and lead to the following physical picture: there are two stable particles in the theory, 
namely the fundamental Schwinger boson fi and the two-boson bound state M2. Higher 
bound states are unstable (resonances) and may decay into all combinations of fi and M2 
final particles that are allowed kinematically. 

There exists also another way for the computation of partial decay widths, namely 
the usual perturbative method of putting all the initial and final states on their mass 
shells, inserting the squared transition matrix element and summing over all kinematically 
possible final states. In our theory this needs the construction of the exact higher n-point 
functions (here the three-point function) as a prerequisite, but once this is done it is easy 
to show that the decay widths (40), (41) just correspond to a first order perturbative 
computation (there the nontrivial numerical factors in (40), (41) are caused by the proper 
normalizations of the initial and final states, which in our computation correspond to the 
residues of the mass poles). 

Observe that the decay widths (40), (41) are related to the M3 binding energy, ~ 
A3. This is a very reasonable result that enables us to interpret the M3 bound state (and, 
a fortiori, higher bound states) as a resonance. Indeed, suppose that the propagator (13) 
belongs to the matrix element of a scattering process (to be discussed elsewhere [§). 
The denominator of (13) has zero real part at s = M| and infinite real part at the real 
production threshold s = = M| + A3, therefore it will give a local maximum of 
the scattering cross section (resonance) near M|, and a local minimum at and the 
resonance width must be related to the binding energy, ~ A3. 
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